Width of the longitudinal magnon in the vicinity of the 0(3) quantum critical point 
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We consider a three-dimensional quantum antiferromagnet in the vicinity of a quantum critical 
point separating the magnetically ordered and the magnetically disordered phases. A specific exam- 
ple is TlCuCls where the quantum phase transition can be driven by hydrostatic pressure and/or by 
external magnetic field. As expected two transverse and one longitudinal magnetic excitation have 
been observed in the pressure driven magnetically ordered phase. According to the experimental 
data, the longitudinal magnon has a substantial width, which has not been understood and has 
remained a puzzle. In the present work, we explain the mechanism for the width, calculate the 
width and relate value of the width with parameters of the Bose condensate of magnons observed 
in the same compound. The method of an effective quantum field theory is employed in the work. 



PACS numbers: 64.70.Tg, 75.10.Jm, 75.40.Gb 



I. INTRODUCTION 

Continuous Quantum Phase Transitions (QPT) and 
behaviour of quantum systems in vicinity of the cor- 
responding quantum critical points attract a great re- 
search interesiji both in theory and in experiment. There 
are manifestations of such effects in cuprate supercon- 
ductors, in iron pnictides and in other materials. QPT 
are most pronounced in low dimensional systems, but 
they also occur in three dimensional (3D) systems. In 
the present paper wc consider quantum systems, there- 
fore 3D means 3D-|-time. A well known material that 
manifests the 3D magnetic quantum critical behaviour 
is TlCuCls. Magnetic and other properties of the com- 
pound has been extensively studied experimentally.—^— 
Under normal conditions the material is nonmagnetic, 
while a magnetic field and/or pressure drives a QPT to 
the magnetically ordered phase. An analogous material 
KCuCla have similar propertiesJ^ii^ It has been under- 
stood that TlCuCla consists of a 3D set of spin dimers 
and the QPT in the magnetic field can be considered as 
Bose condensation of the spin dimer triplet excitations 
(triplons)J.iil"^ The pressure induced QPT can be also 
described in terms of triplonsi^ Inelastic neutron scatter- 
ing from the magnetically ordered phase in the vicinity of 
the pressure induced QPT— clearly indicates two trans- 
verse and one longitudinal spin- wave excitation modes as 
one would expect for condensed triplons. According to 
the datapi^ the longitudinal spin-wave excitation has a fi- 
nite width. Generically this is also expected because of 
the possibility for the longitudinal mode to decay to two 
transverse spin-wave excitations. However to the best of 
our knowledge this width has never been calculated. In 
the present work wc perform such a calculation and com- 
pare the result with the experimental data. The width 
is related to theparameters of the Bose condensate of 
magnons observed in the same compound in the external 
magnetic field. On the technical side in the present work 
we use the effective quantum field theory approach. 

As we pointed out, a description of the system in terms 
of the spin dimer triplon technique is also possible. The 



spin dimer technique accounts for both the short range 
and the long range dynamics on equal footing. On the 
one hand, this is a strength of this description. On the 
other hand, the description is quite technically involved 
if one needs to analyze the interaction effects that we 
are after. The technical complexity is a price for the 
universality. The effective quantum field theory method 
adopted in the present work, allows us to account for the 
interaction effects much more accurately. The method 
does not account for the short range dynamics, but the 
short range effects are irrelevant in the vicinity of a QPT. 

We would like also to note that in the present paper we 
use terms "spin waves" and "magnons" on equal footing, 
having in mind magnetic excitations on both sides of the 
QPT. 



II. EFFECTIVE FIELD THEORY 

The QPT that we consider here is a continuous zero 
temperature transition. Therefore, the short range dy- 
namics are not changed in the QPT, only the long range 
dynamics are changed. In this situation it is quite nat- 
ural to apply the effective field theory approach to de- 
scribe the transitional. A similar approach was used 
to analyse the longitudinal mode in a quasi-one dimen- 
sional antiferromagnet)^ Technically this approach is 
much simpler than the triplon description. The effec- 
tive field theory must be written in terms of the vec- 
tor field 0{t, r) that describes staggered magnetization as 
well as the long-wave-length magnetic excitations. There 
is no net magnetization on both sides of the QPT, only 
the staggered magnetization, hence the effective Landau- 
Ginzburg Lagrangian corresponding to the nonlinear a- 
model describes the system. 
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We note the standard logic for this effective Lagrangian. 
(i) The elastic energy must be proportional to (V:/?)^ 
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(stiffness), (ii) The system consists of spins with the 
same gyromagnetic ratio g = 2. Therefore the Larmor 
theorem must be vahd: If there is a state of the system 
without magnetic field, then the same state in a uniform 
magnetic field B must precess with frequency uj = B 
(Hereafter we set g^B = 1, where is Bohr magne- 
ton.) (iii) The static energy can depend on the magnetic 
field only quadratically because the system has no net 
magnetization. Note that instead of the hard constraint, 
(f^ = const, usually used in the nonlinear cr-model, we 
impose the "soft" quartic interaction ia[(^^]^. The QPT 
is due to the mass term, 

m^ = \'{pc-p). (2) 

In TlCuCla the QPT is driven by external hydrostatic 
pressure p and this is reflected in Eq.(I2|), > is just 
a coefficient. In a more general situation pressure has to 
be replaced by a generalized "coupling constant" . When 
p < Pc and there is no magnetic field, 5 = 0, the mass 
squared is positive and this corresponds to the magneti- 
cally disordered phase with the gapped triple degenerate 
spin-wave dispersion cjg = c^q^ + irfi. When p > pc 
the mass squared is negative and this results in a nonzero 
expectation value of ((/?) that describes the spontaneous 
antiferromagnetic ordering with a gapped longitudinal 
mode and two gaplcss Goldstonc modes. 



III. SPIN- WAVE DISPERSION IN TlCuCla AT 
ZERO PRESSURE 

Since we will compare our results with experimental 
data for TlCuCla we need to determine the correspond- 
ing spin-wave velocity c. TlCuCls has a monoclinic 
structure^ with a unit cell very close to a rectangular 
prism with lattice parameters a = 3.97A, b = 14. 13 A, 
c = 8.87 A and an angle (3 = 96° 10' only slightly devi- 
ating from a right angle (sin/3 = 0.99421...). The spin 
wave dispersion over the Brillouin zone was determined 
by inelastic neutron scattering and the dispersion can be 
very well fit by the following formula— 

efc = V{J + h+9k)^-{fk+9k)^ (3) 

fk = Ja cos kx + Ja2c COs(2fc2: + kz) 
gu = 2Jabc C0s(fc2: -|- fcz/2) cos(fcy/2) 

J = 5.52 mcV 
Ja = -0.24 meV 
Ja2c = -1-57 mcV 
Jabc = 0.46 meV . 

Here — tt < kx,ky < tt, —2tt < k^ < 2tt are di- 
mensionless momenta. The dispersion has minimum 
at fco — (OjO, 27r). Near the minimum the dispersion 
squared can be expanded in powers of momentum, e| = 
-I- aiql + a2ql + a^q'^ + a^qxqz + where q^ = k^, 



qy = ky, qz = kz — 27r, and the coefficients are ex- 
pressed in terms of parameters from Diagonaliza- 
tion of the quadratic form gives the spin-wave velocities 

Cl,C2,C3, 

2 a2, 22, 22, 22 t a\ 

Cfc = ^0 + Cl^l + C292 + C393 (4) 

qi = 0.899g^ + 0.438?^ 

92 = qy 

93 = -0.438g^ + 0.899?^ 
Cl = 7.09 meV 

C2 = 1.12 mcV 
C3 = 0.51 meV 
Aq = 0.65 meV 

The expansion ^ is generally valid only at small mo- 
menta, 91, 92 1 93 ^ 1; however, accidentally the sim- 
ple quadratic dependence along the "softest" direction 
93 is valid practically up to the boundary of the Bril- 
louine zone, 93 ^ tt. Therefore, magnetism in TlCuCls 
is three-dimensional up to the energy /temperature e ^ 
T ~ C37r ~ 1.5 meV fa 15 K above the gap. At a 
higher energy/temperature the magnetism becomes two- 
dimensional. 



IV. SPIN- WAVE EXCITATIONS AT ZERO 
TEMPERATURE AND ZERO MAGNETIC FIELD 

Taking into account Eq.Q the stiffness term ^ ( V(/?)^ 
in the effective action ([Ij has to be replaced by 

\[c\{d,0f + cl{d2^f+cl{d^0f] . (5) 

In the magnetically disordered phase, p < pc the La- 
grangian ([T]) yields the following triple degenerate spec- 
trum 

= + cjql + 4ql + 4ql (6) 

A = Vm^ = X\/pc — P ■ 

Comparing with the gap measured in neutron scattering, 

RefjiS, and presented in Fig. [T] by blue triangles we de- 

1/2 

termine the value of A, A = 0.66 meV/kbar ' . Here we 
take into account that pc = 1.07 kbar. The experimental 
data^ presented in Fig. [T] demonstrate some anisotropy 
that is not accounted in the Lagrangian ([1]). At the mo- 
ment we disregard the anisotropy and return to this point 
later. On the other side of the QPT, p > pc, the field (p 
attains spontaneously the following expectation value 

|<p,| = ./M=^,^-r^. (7) 
V a y/a 

The subscript "c" in tpc stands for classical expectation 
value. We chose the classical field to be directed along the 
z-axis in the spin space. Further standard analysis of the 
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FIG. 1: {Color online) Magnon gaps in the magnetically dis- 
ordered {p < pc = 1.07kbar) and the magnetically ordered 
{p > Pc) phases. The circles, diamonds, triangles etc, show 
experimental data from RefHol. The lines, solid and dashed 
show the theoretical values of the gaps. We take the conven- 
tion that the staggered magnetization in the ordered phase is 
directed along the z-axis, the corresponding longitudinal gap 
is shown by the topmost curve. The Goldstone transverse y- 
magnon is gapless in the ordered phase, while the x-magnon 
is slightly gapped due to the spin-orbit anisotropy. 



effective Lagrangian shows that there are two transverse 
gapless Goldstone x- and y-modes 



= V Cjll + 4<l2 + 4l3 ' (8) 

and one gapped longitudinal z-niode with dispersion 

LJLo - J Al + clql + + 4ql (9) 



Pc 



The "longitudinal" gap agrees pretty well with the data 
shown in Fig. [T] by red triangles and diamonds. 



V. BOSE CONDENSATION IN MAGNETIC 
FIELD AND THE QUARTIC TERM IN THE 
LAGRANGIAN 

To calculate the width of the longitudinal spin wave we 
need to know the coefficient a in the quartic term in (jlj . 
To determine the coefficient we turn to the analysis of 
Bose condensation of magnons in the external magnetic 
field at zero pressure. The Lagrangian ([T|) results in the 
following energy density 

E = 0^4- L 
Sip 

^ l{,n'-B')^' + l{B-0)' + \a[^r .{10) 

Here > since p < Pc- In the energy we skip all 
terms with derivatives since we are interested in Bose 



condensation in a homogeneous system. For magnetic 
field B < m, the energy is minimized when ip = 0, there 
is no spontaneous magnetization. For magnetic field B > 
m (i.e. when the gap is closed), the energy is minimum 
when 



{0-B) = O 

{B^~m^)/a 



(11) 



Thus, there is a spontaneous staggered magnetization 
perpendicular to the direction of the magnetic field. The 
magnetization can be arbitrarily rotated around the di- 
rection of the magnetic field, therefore, the system re- 
tains a U(l) symmetry as it should be for a Bose con- 
densate. Substitution of the expectation value ((TT|) back 
to Eq. pT]]) gives the following energy density of the con- 
densate 



l{B^-r 

E = — 

4 a 



{B — mj 

a 



(12) 



In the second equality in this equation we assume that 
B — m m, i.e. we are close to the phase transition 
point at i3 = TO. This concludes the analysis of Bose 
condensation in terms of the effective action. 

It is well known that the Bose condensation of magnons 
can be also analysed in terms of usual Schrodingcr 
equation within the Hartree-Fock-Popov approximation. 
Such an analysis has been performed earlier in Ref. [22l . 
In the Hartree-Fock-Popov approximation the density en- 
ergy of the condensate at zero temperature reads 



(to — B)nQ + 



vq 2 

— nn 
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(13) 



where no is the condensate density and vq is the effec- 
tive repulsion between magnons. The arrow "— shows 
the result of the energy minimization with respect to uq. 
Comparing (IT2]) and (fT3)) we express the field theory pa- 
rameter a is terms of the Hartree-Fock-Popov parameter 
Vq, a = 2uoTO^. Fit of data4 for Bose condensation in 
TlCuCls performed in Ref. [22| gave the value of the effec- 
tive repulsion, vo ~ 25 meV. The data in Ref. is taken 
at zero pressure, therefore to = Aq ~ 0.65 meV. All in 
all we get 



a = 2woTO^ w 21 meV^ 



(14) 



According to the fit in Ref. [2^ there is a margin of error 
of about 20% in the values of vq and a. 



VI. WIDTH OF THE LONGITUDINAL 
MAGNON 

Decay of the longitudinal magnon is due to the quartic 
term in ([T]) and is described by the diagram shown in 
Fig. m To calculate the decay amplitude one needs first 
to rewrite the Lagrangian in terms of magnons in the 
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FIG. 2: Decay amplitude of the longitudinal magnon. The 
dashed lines indicate magnons and the solid line indicates 
the condensate field tpc- S^pz corresponds to the longitudinal 
magnon, (fi± corresponds to transverse Goldstone magnons. 



magnetically ordered phase. After substitution of (p^ ~ 
ifc + in Eq.([T]) (the condensate field ipc is given by 
Eq.©) we find the quadratic and the cubic part of the 
Lagrangian 



L — L2 + L3 + L4 
1 



(15) 



) C 2 



y(V<5^J^-i|2m2|(5^j2 



L3 = -aifc SLp^{ip±) 

The stiffness term ^ (V...)^ is determined in Eq.([S]). 
Standard quantization of the quadratic Lagrangian gives 
(we set h ~ 1) 



Lk 



"Lk 



'~Px,y 



E 



-Lq 



where a and are corresponding annihilation and cre- 
ation operators. The decay matrix clement shown in 
Fig. [5] is determined after the operators are sub- 
stituted in the cubic Lagrangian L3, 



M = 



(16) 



where q is the momentum of the initial longitudinal 
magnon and qi , q2 are the momenta of the final trans- 
verse magnons. The coefficient 2 in ([TS]) is due to two 
different ways for pairing of creation and annihilation 
operators. The decay width is given by Fermi golden 
rule 
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(27r)3 (27r)3 



X |Mp(5(a;Lq - 



qi 



Lq2)^(9-9l-g2)(17) 



where the coefficient 1/2 stands to avoid the double 
counting of final bosonic states. Integration in PT|) is 



straightforward, the answer is 



a 



87rciC2C3 + c\ql + clql -f- c§g| 

r,=o - ^ Al . (18) 

87rciC2C3 

In this equation (71 , (72 , 93 are the components of the vec- 
tor q. According to the width Fg^o oc cx 
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FIG. 3: [Color online) The decay width of the longitudinal 
magnon at g = 0. The circles and squares, show experimental 
data from RefHol. The solid line shows the theoretical value 
calculated without account of the spin-orbit anisotropy. The 
dashed line shows the theoretical value with account of the 
spin-orbit anisotropy. 



\/p — Pc in a reasonable agreement with experimental 
(lata from RefHol presented in FiglSl According to Eqs. 
(Il8l),(l4]),([l4l) the theoretical value for the ratio F, 

0.21. The corresponding theo- 



g=o/AL 

q=o/^L = s7rciC2C3 ~ The correspondh 

retical curve for F^^q versus pressure is shown in Fig. [3] 
by the black solid line together with experimental datapiS 
shown in the same figure. We stress that we calculate not 
only the functional dependence of the gap on the pres- 
sure. We calculate the absolute value of the gap. Overall 
the agreement between the theory and the experiment is 
quite good having in mind experimental error bars, FiglSl 
as well as 20% theoretical uncertainty in the value of a. 
The agreement can be further improved if one accounts 
for the spin-orbit anisotropy. We discuss the anisotropy 
in the next section. 



VII. ACCOUNT FOR ANISOTROPY INDUCED 
BY THE SPIN-ORBIT INTERACTION 

Data presented in Fig[T] clearly indicate an easy plane 
anisotropy in the magnetically ordered phase: one of the 
transverse magnons is gapped with the constant gap 



0.38 meV 



(19) 
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where the subscript "a" stands for anisotropy. To incor- 
porate the anisotropy in our description wc have to add 
the term 



5Lr, 



1 



2,„2 



(20) 



to the effective Lagrangian ([T]). With account of the 
anisotropy the triple degeneracy of the dispersion in the 
magnetically disordered phase, p < pc, is lifted, the z- 
and y-modcs have the same dispersion (jS]) as before while 
the x-mode has a higher energy 



space. A straightforward calculation shows that Eq. pS)) 
has to be modified in the following way 



rg=o — 



167rciC2C3 



Ar 



(Ai 



4A^) 



(22) 



Here 6{x) is the step function, 9{x) ~ 1 if x > 0, other- 
wise 9{x) = 0. The calculated decay width is shown in 
FigElby the dashed line. The agreement with experiment 
is excellent. 
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The corresponding gaps (wg^o) sue plotted by lines in 
FigUversus pressure. The dispersion in the magnetically 
ordered phase, p > pc, is the following 



X : 



A2 



4<l2 + 4<l3 



y: ujq = J cjqf + 4q^ + clqj 



4<l2 
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A, 



(21) 



The corresponding gaps (w^^o) are also plotted by lines 
in Fig[T] versus pressure. Overall agreement between the 
theory and the experimental data presented in the same 
FiglHis excellent. 

The anisotropy gap influences the width of the longi- 
tudinal magnon only via limitation of the decay phase 



VIII. CONCLUSIONS 

The decay width of the longitudinal magnon in the 
vicinity of the 0(3) quantum critical point of a three- 
dimensional quantum antiferromagnet has been calcu- 
lated using the effective field theory approach. The re- 
sults arc quite generic. 

Wc also compared our results with experimental data 
for TlCuCls where the quantum phase transition can be 
driven by hydrostatic pressure and/or by external mag- 
netic field. The parameters of the effective action were 
determined from available data for inelastic neutron scat- 
tering and for Bose condensation of magnons in the exter- 
nal magnetic field. Having determined these parameters, 
we calculated the decay width. The result agrees very 
well with measurements. 
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